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ABSTRACT 

Context. After the discovery of the first connection between GRBs and SNe almost two decades ago, tens of SN-like rebrightenings 
have been discovered and about seven solid associations have been spectroscopically confirmed to date. 

Aims. The luminosity and evolution origin of three SN rebrightenings in GRB afterglow light curves at z ~ 0.5 are determined 
along with accurate determinations of the host-galaxy extinction. Physical parameters of the SN explosions are estimated, such as 
synthesised 56 Ni mass, ejecta mass, and kinetic energy. 

Methods. We employ GROND optical/NIR data and Swift X-ray/UV data to estimate the host-galaxy extinction by modelling the 
afterglow SED, to determine the SN luminosity and evolution, and to construct quasi-bolometric light curves. The latter were corrected 
for the contribution of the NIR bands using data available in the literature and blackbody fits. Arnett's analytic approach has been 
employed to obtain the physical parameters of the explosion. 

Results. The SNe 2008hw, 2009nz, and 2010ma observed by GROND exhibit 0.80, 1.15, and 1.78 times the optical (/' band) 
luminosity of SN 1998bw, respectively. While SN 2009nz exhibits an evolution similar to SN 1998bw, SNe 2008hw and 2010ma 
show earlier peak times. The quasi-bolometric light curves (340-2200 nm) confirm the large luminosity of SN 2010ma (1.4 x 10 43 
erg s -1 ), while SNe 2008hw and 2009nz reached a peak luminosity closer to SN 1998bw. The modelling resulted in 56 Ni masses of 
around 0.4-0.5 M 0 . 

Conclusions. By means of the a very comprehensive data set, we found that the luminosity and the 56 Ni mass of SNe 2008hw, 
2009nz, and 2010ma resembles those of other known GRB-associated SNe. This findings strengthens previous claims of GRB-SNe 
being brighter than type-Ic SNe unaccompanied by GRBs. 

Key words, gamma-ray burst: individual: GRB 081007 - supernovae: individual: SN 2008hw - gamma-ray burst: individual: GRB 
091127 - supernovae: individual: SN 2009nz - gamma-ray burst: individual: GRB 101219B - supernovae: individual: SN 201 Oma 


1. Introduction 


Gamma-ray bursts (GRBs) and supernovae (SNe) correspond to 
the most energetic explosions in the Universe with a radiative 
energy release of about 10 51-53 erg. Nowadays the observational 
evidence points towards the catastrophic deaths of massive stars, 
w hich are thought to give birth to both long GRBs (durations ^ 2 
s: iKouveliotou et all 19931 ) and broad-lined (BL) type-Ic SNe af- 
ter the collapse of their cores into a black hole ( BH; iPacz vnski 
1 998aL [Fryer et alii 9991: van Paradiis et al.ll2000l). Known as the 


collapsar model (Wooslevl|l993b iMacFadven & Wooslevl[l999l : 

Bromberg et all2012 ), the collapsing core of a very massive star 
can lead to the formation of a relativi s tic jet that will produce 
high-energy emission dWooslevi [ 19931 : iWooslev & MacFadvenl 
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1999) in the form of a GR B or an X-ray flash (XR F; lHeise et al.l 
2001 : iKippen et all 12004 ISakamoto et al.1120081) . The y-ray 
emission itself lasts from a few tenths of a second to a few thou¬ 
sand seconds, is generated within the outflow at ultra-relativistic 
velocities, and is collimated into a jet dZhang et alJ 12009 ) that 
drills its way out of the star. The interactions between fireball 
shells with different speeds (“internal shocks”) are responsible 
for the prompt y -ray emission. The multi-wavelength afterglow 
(AG), detectable fro m radio throughou t to X- rays up to months 
after the GRB (e.g., iKann et aT1l2010t 1201 ll) . is explained by 
the synchrotron emission produced in the interaction between 
the circumb urst material and the rela tivistic outflow (“external 
shocks”; see lZhang & Meszarosil2004l for a review). 

In principle, the energy transferred to the envelope should 
also be capable of causing the ej ection of the stellar envelope 
dBurrowsll2000t iHeger et aDl2003l) . However, it is unclear how 
or even if there is always enough energy for the SN explo¬ 
sion (“fail-back” events, e.g., iFrver et alll2007i . and references 
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therein). Moreover, it is unknown exactly to what extent the pro¬ 
genitors have to lose their envelope to produce a GRB. However, 
it is generally accepted that type-lb and type-Ic SNe are formed 
from evolved high-mass progenitors like Wolf-Rayet (WR) stars, 
which have liberated their outer shells through (1) pre-SN stellar 
winds, (2) mass transfer to a binary companion due to Roche- 
lobe overflow, or (3) a combination of both processes. The stel¬ 
lar explosion is then referred to as a “stripped-envelope” SN (SE 
SN; lClocchiatti & Wheelei][T99% . The end result of a GRB-SN 
explosion wou ld correspond to a comp act remnan t, either a neu - 
tron star (NS; iBaade & Zwickvl 119341) or a BH (I Arnetll 1 1 9961) . 
To date, long GRBs have been associated only with type-Ic BL 
SNe, which are th ose lacking H (Minkowsk il 19411) and He lines 
( Filinnen kol 19971) and showing expansion velocities in the order 
of 20,000 km s -1 (for re vie ws on the GRB-SN con nection, see 
IWooslev & Bloonill20061 and iHiorth & Bloomll2012l) . 


The first and most representative case of the GRB-SN 
connection is the associ ation of SN 1 998bw with the u nder- 
luminous GRB 980425 (lKippenlll998l: ISadler et al.lll998l). Al¬ 
thoug h initially controversial ( Galama et al.1 19981: Pian et all 
1998), the physical association between t hese events was sup- 

and spatial grounds 
. Five years later. 


ported on t empor al dlwamotojhaL 


dPian et akl 1200(1 iKouveliot ou et al 


99 


200 . 


the association of GRB 030329 with SN 2003dh was clearly 
identifie d through spectra showing both the AG and SN coun¬ 


terpart s (IHiorth et al 


120031: iMatheson et al 


2003t iKawabata et~aTl l2003t IStanek et~aT1 
120031) and became a solid piece of evi¬ 


dence in favour of the GRB-SN connection. There have been 
a number of other spectroscopic associations]]], which in the lit¬ 
erature are_also dubbed “hypernovae” (HNe; Paczvn ski 1998b; 
lHansenlfl999l) given their high luminosities. Due to their high 
energetics, HNe produce ^ 0.2 M 0 of 56 Ni, are thought to have 
very massive progenitors , and are often connecte d to BH forma¬ 
tion dNomoto et al.ll2010l:IStritzinger et al.ll2009l) . 


Late-time rebrightenings in AG light cu rves have been inter- 


preted as SN signals, e.g., GRBs 970228 ( 

Galama et al.| 2000|: 

Reichartet al. 2000), 011121 dBloometal. 

2002; Greiner et al. 

2003), 020405 (Price et al. 2003; Masetti et al. 2003), 041006 

(Stanek et al.l 2005 [ Soderberg et al. 

2006b, 060729, and 090618 

(the latter two in Cano et al.l 20 ll 

) to mention a few. These 


photometric bumps are consistent in terms of colour, timing, 
and brightness with those expected for the GRB-SN popula¬ 
tion (IZeh et al.ll2004t iFerrero et al.ll2006l) . but they are usually 
at faint apparent magnitudes, which hampers the spectroscopic 
identification. However, the SN counterpart can be as b right as 
My = -19.8 mag for SN 20031w dMalesani et al.l2004 ). These 
rebrightenings have been detected in AG light curves out to red- 
shifts of ~ 1 (e.g.. iMasetti et al] 120051: [Delia Valle et al.l [20031 
owing to the sensitivity of current ground-based telescopes ded¬ 
icated to follow-up observations. A handful of sample studies 


1 These also showed BL features in their spectra: GRB 021211/SN 
20021t dDella Valle et all 120031 ). GRB 031203/SN 20031w (M alesani 
et al. 2004), GRB 050525A/SN 2005 n c dDella Valle et alj|20 06). GRB 


060218/SN 2006ai dPian et all 120061: iModiaz et all l200fil: Sollerman 


et al. 2006), GRB 081007/SN2008hw {Del la Valle et al.||2008| : Jin et al. 
2013), GRB 091127/S N 2009nz dBereer et alfeoil ). GRB 

ID, 


101219B/SN 2010ma ( Sparr e~etal] |201 il l. GR B 1 I 121 I Aide 


Postigo et al 2012 ), GRB 100316D/SN 2010bh dChomock et al 


Ugarte 


2010 ; 


Bufan o et alJUOlZ) GRB 120422A/SN 2012bz (Melandri et al.l 12012 ; 
Schul ze et all l2014l).GRB 120714B/SN 2012eb klose et all2012allbl) . 
GRB 130215A IdTusarte Postigo et al 112013l:ICano et alJl2014l) . GRB 
130427A/SN 2013c q dXu et all l2013li . GRB 13 0702A/SN 2013dx 
dSchulze et all l20l3l) . GRB 130831 A/SN 2013fu dKlose et all 11)11 
ICano et all2014l) . and GRB 140606B dPerlev et al.ll2014l) 


of GRB-SNe (including bumps not spectroscopically identified) 
have analysed the luminosity distribution, the light-curve mor¬ 
phology, and the explosion physical parameters such a_s_ kinetic 
j y (Fk), ejected mass ( M P ,), and 56 N i mass ( Mm; : iRichardsonl 
iThiine et~akl 1201 It ICanol 120131) . They concluded that 
GRB-SNe are in general brighter than the local sample of SE 
SNe, except for cases of exceptio nally bright ty pe-Ic SNe (e.g., 
SN 2010av. [Sanders et al.l 1201 2t SN 201 0gx, [Pastorel lo et al.l 
1201 Oh . Regarding l ight-curve mo r pholog y, IStanek et al.l (2005) 
and more recently ISchulze et al.l ( 2014 ) claim to have found 
a correlation between brightness and light-curve shape, which 
was also confirmed by ICand (12014l) using a larger sample and 
including the appropriate K corrections. This strengthens the 
use of GRB-SNe as standard candles for cosmology (see also 
recent studi es bv iLi & Hiorthl 120141 ICano & Jakobssonll2014l 
and ILi et al.l 120141) . While more than two dozen photometric 
bumps in AG light curves have been presented as SN rebright¬ 
enings (e.g.. lRichardsonl2009l) . so far only seven have been con¬ 
firmed through high-S/N spectra: SNe 1998bw, 2003dh, 20031w, 
2006aj, 2010bh, 2012bz, 2013dx, and 2013cq. 

The energy injection of a newly-formed rapidly-spinning 
strongly-magnetised NS (so-called “magnetar”) provides an al¬ 
ternative scenario for GRB-SNe. Here the SN is powered 
by the dipole-field stren gth of the magnetar (e.g., IWooslevI 
1201 Ot: iDessart et al.ll2012h . Magnetars have been linked to the 
GRB emission too, because their outfl ows can explain the en¬ 
ergetics of long-du ration GRBs (e.g., iBucciantini et al.l [2009; 
iMetzger et alJl20lTh . Moreover, the E k of GRB-SNe (~ 10 52 
erg) is fairly consistent with the maximum rotational energy of a 
NS with a period of 1 ms (iMazzali et al.l2006ah . GRB-SN zoo is 
claimed to be entirely produced b y mag netars and driven by the 
SN rather than by the GRB jet (IMazzali et al.ll2014l) . 

Three detected SNe associated with GRB counterparts are 
the main focus of this paper: SNe 2008hw (GRB 081007), 
2009nz (GRB 091127), and 2010ma (GRB 101219B). The ac¬ 
quisition, reduction, and calibration of the multiwavelength data 
are described in Sect. [2] The corresponding analysis is presented 
in Sect. [3] along with further discussion in Sect. [4] Finally, we 
summarise our conclusions in Sect. 0 


2. Data 


For the three objects of interest the data was obtained by the X- 
Ray Telescope (XRT; iBurrows et al. 20 05h and the Ultra-Violet 
Optical Telescope (UVOT; Romin g et al.l 120051) both on board 
the Swift satellite (iGehrels et al.l 2004 ) and by the Gamma- 
Ray burst Optica l and Near-infrared Detector (GROND; 
iGreiner et al.l 120071 120081) . the seven-channel imager mounted 
on the MPG 2.2-m telescope at La Silla, Chile. The whole data 
set comprises X-ray photometry and spectra from 0.2-10 keV, 
UV/optical photometry in the uvw2 uvml uvwl ubu filters, and 
optical/NIR photometry in the g' r' i' z' J H K s bands, spanning 
four orders of magnitude in the energy spectrum. 

The UVOT/XRT data retrieval and the GROND/UVOT 


methodology towards the final photometry are detailed in 
lOlivares E. et al.l ( 2012 j). Optical image subtraction of the host 
galaxy was performed for GRB 081007/SN 2008hw and GRB 
091127/SN 2009nz. All data presented are corrected for the 
Galactic foreground reddening E(B - V)Gai in the direction of 
the burst (ISchlegel et al.l 119981) . The reddening is transformed 
to the extinction Aycai by assuming a ratio of total to selective 
absorption of R i:c, a i = 3.1 from the Milky-Way (MW) reddening 
law. The final GROND photometry is presented in Appendix [Al 
All magnitudes throughout the paper are in the AB system. 
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Table 1. GROND sample of GRB-associated SNe. 


GRB 

SN 

RA(J2000) 

|- h . m . s j 

Dec.(J2000) 

[ °: ': " ] 

z 

d & 

[Mpc] 

^V,Gal b 

[mag] 

NH,Gal C 
[10 2() cm” 2 ] 

Refs. 

081007 

2008hw 

22:39:50.40 

-40:08:48.8 

0.530 

2885 

0.05 

1.4 

1 

091127 

2009nz 

02:26:19.87 

-18:57:08.6 

0.490 

2628 

0.12 

2.8 

2 

101219B 

2010ma 

00:48:55.35 

-34:33:59.3 

0.552 

3022 

0.06 

3.1 

3 


Notes. ( ^_The luminosity distances are computed using the ACDM cosmolog ical m odel (O m = 0.27, fi A = 0.73, and H 0 = 74.2 km s~* Mpc~'; 
iRiess et al.l2009h and the redshifts corrected by the local velocity field (Mould et al. 2000). <b) The Galactic foreground e xtinction values are ta ken 
from the dust maps of Schlegel et al. (1998). (c) The absorption column densities are taken from the Galactic Hi maps of Kalberla et al. (20051). 

References. Redshifts are taken from ( ll lBereer et al .1120081 . 12l lVergani et al.ll201 ll Gi lSnarre et al.ll2011bl . 

Table 2. GROND photometry of the host galaxies. 


GRB 

SN 

9' 

r' 

V 

z' 

J 

H 

K s 

081007 

2008hw 

24.66 + 0.11 

24.49 + 0.11 

24.08 + 0.19 

23.96 ± 0.24 

> 22.0 

>21.1 

>20.1 

091127 

2009nz 

24.08 ± 0.09 

23.45 + 0.06 

22.85 + 0.07 

22.95 ± 0.08 

> 21.7 

>21.4 

> 19.9 

101219B 

2010ma 

> 25.4 

> 25.2 

> 24.5 

>24.5 

> 22.2 

> 22.0 

> 20.2 


Notes. The host-galaxy magnitudes are all corrected for the corresponding Galactic foreground extinction. The upper limits were derived from the 
deepest observation available showing no detection and are quoted at the 3<x confidence level. 


3. Three GRB-associated SNe detected by GROND 

TableQ]presents a sub-sample of GRBs with late-time optical SN 
rebrightenings in their AG light curves, all of them observed by 
GROND. Deep late-time observations were carried out for each 
of them to constrain the contribution from their host galaxies. 
If the host was detected, we performed image subtraction. Ta¬ 
ble |2] presents the resulting photometry for those host galaxies. 
In the following, observational facts and general properties of 
each event are summarised from the literature. If possible, mass 
estimates are deri ved from the SEDs o f the host galaxies using 
the hyperZ code (fBolzonella et al.l[2000l) and a library of galaxy 
spectral templates extinguished by the different reddening laws. 


GRB 081007/SN 2008hw The Swift /BAT (iBarthelmv et al.1 
2005) discovered GRB 081007 at 05:23:52 UT on 2008 Octo¬ 
ber 7 ([Baumgartner et alll2008i) . The prompt emission had a du¬ 
ration of 71) 0 ~ 10 s and a soft spectrum with /: pca k S) 30 keV 
(iMatTwatr^eDdl 200 I). The redshift of z — 0.5295 was found 
bv iBerger et al] ( 200 8) through optical spectroscopy. A subse¬ 
quent optical spectrum taken 17 days after the burst shows broad 
features indicative of an emerging SN, which was thereafter clas¬ 
sified as type I (no Hydrogen lines) and named SN 2008hw 
dDella Valle et al.l[2008l) . The SN bump was also reporte d as a 
flux excess with respect to the AG (Soderberg et al. 2008). The 
GROND photometry of the host galaxy (Table [2} from August 
31, September 30, and October 21, 2011, yields a stellar-mass 
range of M* ~ 10 8 ~ 9 M 0 , which is compatible with the popu¬ 
lation of GRB hosts dSavaglio et al.l 1200 9). Using appropriate 
transformation equation^ our host magnitudes are somewhat 
brighter but marginally consistent with the measurements pub¬ 
lished bv lJin et akl d2013l) of Rc > 24.67 and Ic = 24.29 ± 0.20 
mag at ~ 87 d after the GRB. 


GRB 091127/SN 2009nz At 23:25:45 UT on 2009 November 
27, the Swift/BAI was triggered by GRB 091127 (ITroia et al.l 
l2009h . The y-ray emission lasted for Tgo = 7.1 s and 
showed a soft spectrum (IStamatikos et al.l 12009: Troia et all 


[20l2h . A redshift of z_ = 0.490 was obtained fr om op - 
tical spectroscopy (ICucchiara et al.l l2009t iThone et~aD [2009). 
Observatio ns by Konus-Wind con firmed the results from the 
Swift/BAI dGolenetskii et al.l|2009l) and additionally yielded an 
energy release typical for cosmological GRBs (£ r ,i SO ~ 10 52 
erg). The optical AG was confirmed with GROND observations 
(lllpdike et al.l l2009) adding NIR detections. The full ana lysis o f 
the GROND AG light curve was presented in lFilgas et al.l(l201 ll) . 
The SN classification became official based on the photomet¬ 
ric SN bump (Cobb et al. 2010a bj) and spectroscop y was pub¬ 
lished later showing BL features dBerger et al.ll201 ll). Photome¬ 
try dep icting the SN rebrightening was published in lCobb et alJ 
(201 Og) and IVergani et all d201 ll) . Using the host-galaxy de¬ 
tections in GROND optical i mag ing (October 31, 2010) and in 
NIR photometry from IVergani et al.l d201 ll) . a stellar mass of 
M* = 1O 84 M 0 is obtained. This value falls in the low- mass end 
of the observed distribution of GRB host masses dSavagflo et al.l 
2009) and is comp atible with the stellar mass computed by 


Vergani et al.l d201 ll) . 


http://www.sdss.org/dr7/algorithms/sdssUBVRITransform 


GRB 101219B/SN 2010ma At 16:27:53 UT on 2010 Decem¬ 
ber 1 9, the Swift/BAI discovered GRB 1 01219B dGelbor detalJ 
l2010t) . The BAT burst lasted Tgo ^ 34 s ([Cummings et al,ll2010t) 
and consisted of a spe c trum with £p e ak — 70 keV as observed 
by Fermi /GBM dvander Horst! 1201(1). The SN d i scove ry was 
first reported photometrically bv lOlivares E. et akl (1201 ll) along 
with a redshift estimation assuming the brightness of SN 1998bw 
for the rebrightening (z. = 0.4-0.7). The spectroscopic confir¬ 
mation of SN 2010ma came later bv ISparre et al.l d201 ll) along 
with the redshift determination of z - 0.55185 from weak Mg 
absorption lines. The spectroscopy lead to further analysis by 
ISparre et al.l (1201 lbt) that shows broad-line features characteris¬ 
tic of GRB-SNe. Late-time GROND observations on September 
30, 2011, show no signal of a host-galaxy down to deep lim¬ 
its (Table O, therefore no image-subtraction procedure was per¬ 
formed. These upper limits imply a stellar mass for the host 
galaxy of M* < 1 0 9 ' 2 M Q , which corresponds to the low-mass 
half of observed GRB host mass distribution and is marginally 
bomrhpatible with the Small Magellanic Cloud (SMC). 
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3.1. Multicolour light-curve fitting 


After image subtraction of the host galaxy in the cases where 
it was detected (Table [2}, the light curves were fitted simul¬ 
taneously using one or two power-law components (F v ~ t a ) 
and templates of SN 1998bw, where corrections due to redshift 
and Galactic foreground extinction were taken into account (see 
IZeh et all 12004 for details on the fitting of SN 1998bw tem¬ 
plates). Simultaneous modelling consists of unique power-law 
slopes a i, 0 - 2 , and SN-template stretch factor .v for all bands. 
The ratio between the luminosity of the observed SN and that of 
SN 1998bw (luminosity ratio k: IZeh et~aT]|2004:) represents the 
free brightness parameter, which was fitted to the light curves 
corrected for Galactic extinction only. Therefore, the luminosity 
ratios are then corrected for the host-galaxy extinction A y^ost de¬ 
termined by the SED modelling (see Sect. 13. 2I >. The modelling 
is described in detail below for each event and summarised in 
Tables [3] and [4] 


GRB 081007/SN 2008hw Figure Q] shows that the light curves 
in all seven bands are w ell modelled using a broken power law 
dBeuermannet al.1 fl999). The g'r'i'z' photometry has been 
image-subtracted to remove the host-galaxy flux. The X-ray 
light curve from the Swift[X. RT was included in the fitting to 
constrain the decay after the break, where there is only a sin¬ 
gle optical epoch. For the r'i'z' bands, it was necessary to add 
a supernova component with a luminosity about 65 - 80% that 
of SN 1998bw (see Table Q}. Due to the JHK s flux excesses 
with respect to the broken power law at roughly 1 day after the 
burst, a constant component was included in the modelling for 
these bands. The g -band upper limit is strongly affected by 
absorption of metal lines and wavelength extrapolation of the 
SN 1998bw template (e.g., the case of SN 2009nz due to high 
redshift). l.Tin et al. i! ( 2013 ) report a luminosity 50% that of SN 
1998bw, however, without accounting for the significant host ex¬ 
tinction (see Sect. 13.2k 


GRB 091127/SN 2009nz Figure [2]presents the light curves of 
the AG in six bands. The g'r'i'z' photometry has been image- 
subtracted to remove the host-galaxy flux. All are well fitted 
by a single power law, which needed a SN component for the 
g'r'i'z’ bands. No X s -band observations were obtained for this 
event (iFilgas et al.ll201 ll) . The brightest host galaxy allowed by 
the data was included in the model for JH at late times (see Table 
0. The k and .? values reflect strong similarities to SN 1998bw in 
the r'i 'bands. At the redshift of SN 2009nz (z — 0.490), the g' 
band is probing wavelengths centred at ~ 3000 A, where the flux 
is strongly affected by absorption-line blanketing of metals, and 
so the intensity can differ from SN to SN. Moreover, since the U 
band, the bluest band from which the SN 1998bw templates are 
constructed, is sensitive ^ 3000 A, extrapolations dominate the 
g -band template. Given also the non-detections after day 12, we 
derived an upper limit of k g > < 1.21 from the fitting (Table [4j». 

GRB 101219B/SN 2010ma Figure 0 shows the GROND light 
curves of the optical counterpart. The SN bump is clearly seen 
in the r’i'z’ bands, however, it is less significant in the g' band. 
At the redshift of the event, the gr'band actually probes the UV 
regime, therefore, the lower g -band SN luminosity is explained 
by a combination of both the wavelength extrapolation of the 
template and the UV line blanketing by metals. Even though 
the host galaxy remained undetected, it may explain the flux ex- 
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Fig. 1. Multicolour light curves of GRB 081007/SN 2008hw corrected 
for Galactic extinction as observed by the SwiftPiKT (upper panel) and 
GROND (lower panel). Filled circles represent detections and arrows 
are upper limits. Solid lines correspond to the overall fits and dotted 
lines to the AG component. For clarity, light curves were shifted along 
the magnitude axis. Shallow upper limits are not shown (see Table I A. fl 
for the complete data set). 



Fig. 2. Multicolour GROND light curves of GRB 091127/SN 2009nz 
corrected for Galactic extinction. Only data after day one has been em¬ 
ployed. The symbol and line coding is the same as Fig.[7]as well as the 
vertical shift for clarity. 
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Table 3. Parameters of the AG component and goodness of the light-curve modelling. 


GRB 

at 

f break [days] 

1 

1+2 

X 2 /P 

081007 

-0.66 ± 0.01 

0.91+0.05 

15 fixed 

-1.40 + 0.05 

1.5 

091127 

-0.38 + 0.01 3 

0.34 + 0.0P 

1.3 + 0.1 3 

-1.63 + 0.02 

1.4 

101219B 

-1.01+0.01 




1.8 


Notes. The primary power-law slope is at . In case of a break in the light curve, a secondary slope 02 along with the break time fb rea k and break 
smoothness parameter 77 are introduced. The ratio x^/fJ. * s computed in the multiple-component fitting procedure, which includes AG plus S N 
modelling. See Table[4]for the SN parameters. (a) Parameters were taken from the fitting of the full GROND r -band light curve bv lFilgas et alJ 
d201 ll) except a 2 , which was fitted by a single power law using the data presented in Fig.[2]only. 

Table 4. Parameters of the SN component with respect to SN 1998bw templates. 


SN 

GRB 

Stretch 
factor (s) 

9' 

Luminosity ratio (k f 
r' i' 

z' 

2008hw 

2009nz 

2010ma b 

081007 

091127 

101219B 

0.85+0.11 

1.03+0.04 

0.76 + 0.10 

<0.90 
< 1.21 

0.85 + 0.17 

0.80 + 0.10 
1.15+0.09 

1 70 +0.08 
^ -0.17 

0.65 ± 0.08 

0.96 + 0.14 

1.36 + 0.09 

0.69 + 0.10 
0.73 + 0.12 
0.63 ± 0.09 


Notes. <a) Luminosity ratios are all corrected for Galactic and host-galaxy extinction. The latter correction is taken from the AG SED fitting in 
Sect. 13.21 <b) No host-galaxy contribution was assumed. See text for estimations including host emission. 
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Fig. 3. Multicolour GROND light curves of GRB 101219B/SN 
2010ma corrected for Galactic extinction. The symbol and line cod¬ 
ing is the same as Fig.[T|as well as the vertical shift for clarity. The red 
dashed line represents a model with an extra host-galaxy component. 


cess 35 days after the burst in the r' band (dashed line in Fig. 
|3}- The k value would decrease ~ 14% in this case. Therefore, 
the lower error in k r ' was increased to match the 3 cr lower limit 
when assuming the brightest host component possible (Table |4j. 
We note that the k value for the z' band is smaller compared to 
the bluer bands. Along with the differences in the SN luminos¬ 
ity ratio among all bands, this indicates that the colours of SN 
2010ma are different from those of SN 1998bw. 

Fig. [4] shows the colour curves of the three SNe analysed 
compared against the templates of SN 1998bw, where the bluer 
emission of SN 2010ma is significant at early times. 


Fig. 4. Colour curves corrected for the total extinction of SNe 2008hw 
(blue circles), 2009nz (purple squares), and 2010ma (gold diamonds) 
after AG and host subtraction. Blue, purple, and gold solid lines are 
computed from the templates of SN 1998bw at the redshifts of SNe 
2008hw, 2009nz, and 2010ma, respectively. 


3.2. Spectral energy distributions 

Using the available X-ray data from the Swift/X RT, the 
UV/optical data from the SwiftfUV OT, and the optical/NIR data 
from GROND, we constructed a single AG SED per event with 
the main purpose of determining the extinction along the line 
of sight through the host galaxy. The SED modelling was per¬ 
formed similarly as in [Greiner et al.l d201 If) and the results are 
presented in Table 0 
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Table 5. Parameters of the SED modelling of the AG. 


GRB 


A l/host 

[mag] 

^break 

[eV] 

-VH.host 

[10 21 cnT 2 ] 

x 2 /v 

081007 

0.97 ± 0.09 

0.68 ± 0.08 (SMC) 

27 + 54 

J1 -12 

5.6 ± 0.7 

1.1 

091127 a 

0.748 ± 0.004 

< 0.03 (LMC) 

2.6-29.9 

0.32 ± 0.06 

1.1 

101219B 

1.12 + 0.01 

0.12 + 0.01 (SMC) 

9.0 fixed 

0.6 ± 0.3 

0.8 


Notes. Obeying the fireball model for GRB AGs, the high-energy (fix) and the low-energy (0 opt ) spectral indexes are correlated by fix = /3 0p t + 0.5 
for Ebrcak ~ ^cooling, where the latter comes from the cooling frequency of the electrons, except in the case of GRB 091127 were /f opt varies in the 
range 0.25^0.62A a) The quoted values of fix, (V H ,host. and reduced^ 2 are computed from the simultaneous best fit to all ei ght G ROND/XRT SEDs 
bv iFileas et alJ 1 20111 1. The Ebreak range comes from an observed evolution of /3 opt . The A y^ost upper limit was taken from lSchadv et all ( 20120 . 


Note that the AG may probe a slightly different line of sight 
than the SN photosphere. If anything, the extinction for the SN 
should be larger than for the AG, because the AG forms further 
out, where the material ejected by the GRB hits the circumstel- 
lar medium. In the standard fireball shock model, this radius 
is ab out 10 17 cm (and even larger for low-luminosity events; 
iMolinari et al.ll2007tl . Moreover, dust can be formed in the SN 
ejecta, although not significant amounts on such short timescales 
(e.g.. [Smith et aljfeOlA . Therefore, we considered that the ex¬ 
tinction determined through the AG SED is valid for the SN 
component as well. The following corresponds to a description 
of the SED fitting for each of the rebrightenings. 


GRB 081007/SN 2008hw To include contemporaneous 
SwiftfUVOT data, the second GROND epoch was chosen to 
study the broad-band SED of GRB 081007 presented in Fig. 
[3 From the UVOT, upper limits in the UV bands are included, 
which help to constrain the host-galaxy extinction. The time- 
integrated Swift/X RT spectrum was interpolated to the epoch 
of the UV/optical observations. The resulting values of host- 
galaxy extinction and their correspondin g stat istical uncertainty 
are consistent with those computed by ICovino et al.l d2013l ) and 
for the GROND filters we obtain A g -host = 1.39 ± 0.16, A- '.host = 
0.99 ± 0.12, A,-host = 0.77 ± 0.09, A-'^ost = 0.63 + 0.07, Ayhost = 
0.38 +0.04, A/yhost = 0.24 + 0.03, and A^.host = 0.14 ±0.02 in the 

observer’s frame, all in units of magnitude. These values were 
used to correct the SN luminosity ratios shown in Table [4] 


GRB 091127/SN 2009nz The broad-band SEDs of the early 
AG of GRB 091127 were presented bv lFilgas et al. ( 201 ll) using 
the GROND data. A detailed analysis by iSchadv et alJ ( 2012 1 
includes Swift/(JV OT and GROND data and constrains the host- 
galaxy extinction, which results in Ayhost < 0.03 mag. The SED 
parameters are shown in Tabled 


GRB 101219B/SN 2010ma Using GROND, XRT, and UVOT 
data combined, the AG SED of GRB 101219B was constructed 
at 9 h after the burst. Figure [6] shows a broken power law as the 
best fit. The values of the required host-galaxy extinction for the 
GROND filters and their corresponding statistical uncertainty in 
the observer’s frame are A^host = 0.25 ± 0.03, A,',h os t = 0.18 ± 
0.02, A, 'j 10 st = 0.14 ± 0.02, A ; 'host — 0.11 ± 0.01, A yhost = 0.07 ± 
0.01, A ft host = 0.04 ± 0.01, and A^host = 0.026 ± 0.003, all in 
units of magnitude. 



Energy [keV] 


Fig. 5. Broad-band AG SED of GRB 081007 at 1.6 ks after trigger. 
The arrows are 3cr upper limits. The best-fit model (thick line) is an 
extinguished broken power law. The thin line represents the unextin¬ 
guished model. The residuals are in units of^ (lower panel). 



Energy [keV] 


Fig. 6 . Broad-band AG SED of GRB 101219B at 9.0 h after trigger. 
The symbols, the line coding, and the panels are the same as in Fig. [3 
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3.3. Calculation and modelling of the bolometric light curves 

To isolate the SN from the AG evolution, the light-curve mod¬ 
els computed in Sect. 13.11 were employed. The AG contribu¬ 
tion was calculated from the model for the epochs when the SN 
bump was observed and it was subtracted from the light curves 
for each filter. The uncertainties in the model were appropriately 
propagated to the final magnitude errors. After the AG subtrac¬ 
tion, quasi-bolometric light curves were computed for each of 
the three events by numerically integrating the monochromatic 
fluxes in the wavelength range from 340 to 700 nm. The redshift- 
based luminosity distances in Table [T| were employed to trans¬ 
form observed into absolute flux. The total uncertainty in the 
luminosity distance is about 10% and has not been included in 
the quasi-bolometric light curves. 


3.3.1. NIR bolometric correction 


The NIR luminosity proves critical when estimating the bolo¬ 
metric flux and consequently the physical parameters of the ex¬ 
plosion obtained via the quasi-bolometric flux. However, SNe 
2008hw, 2009nz, and 2010ma remained undetected in the JHK 
bands. To account for the NIR flux of these SNe, we proceeded 
with two different methods. First of all, we defined the NIR flux 
from 700 to 2200 nm in the rest frame and the quasi-bolometric 
flux from 340 to 2200 nm. The quantity to be estimated via the 
two methods is the ratio between the NIR flux and the quasi- 
bolometric flux as defined above. 

The first method consisted in estimating the NIR fraction of 
the quasi-bolometric flux using the observed NIR data available 
in the literature. With the optical/NIR photometry, we com¬ 
puted the ratio between the NIR and the quas i-bolometric fluxes 
for the GRB-SN e 1998bw (iKocevski et alJ l2007h and 2006aj 
jPatatemdj200]k and for the type- Ib/c SNe 2002ap dFolev et al'l 
120031: lYoshii eHfl1f20()l . 2007uy dRov et al.ll20i3h . and 2008D 
(iModiaz et al.l2009l) . We show the observed NIR fractions in the 
upper panel of Fig. [7] along with a quadratic-polynomial fit for 
each SN. For each fit, we also obtained the corresponding uncer¬ 
tainty as a function of time. By taking the weighted average per 
time bin (5-days width, f ( - = 1 d), we derived the joint evolution 
of the NIR fraction for the five SNe. The non-weighted RMS 
was taken as the lcr error. Only the first time bin uses data from 
a single event and the error here is approximated by the uncer¬ 
tainty in the individual polynomial fit. The binned NIR fraction 
was interpolated using quadratic polynomials to retrieve values 
for a given time and to plot the grey contours in Fig. [7] The SNe 
2007uy and 2008D do not contribute much to the weighted av¬ 
erage, because their host extinction is large an d highly uncertain 
(0.3-0.5 mag: lMazzali et al.ll2008HRov et al.ll2013h . 

The second method assumes that the SN atmosphere at early 
phases resembles a cooling black body ( BB; e.g.. I ArnetBl 1 982h 
FiliD nenkolll997l:rDessart & Hillierll2005r> . We defined a simple 
BB model with the temperature and a flux normalisation as free 
parameters. Then, we modelled the SN SEDs constructed using 
g'r'i'z' data at each epoch. The results of the fitting proce¬ 
dure are shown in Appendix O We obtain colour temperatures 
decreasing with time and consistent with other SE SNe (e.g., 
iFolatelli et al.ll2014 i). The extrapolation into the NIR range de¬ 
livered NIR fractions plotted in the lower panel of Fig. [7] with 
uncertainties between 0.05 and 0.13. The increasing NIR flux 
with time is consistent with the scenario of the cooling envelope. 
We repeat ed this procedure fo r the optical data of SNe 1998bw 
(UBVRI ; iGalama et alJ 19981) and 2006aj ( UBVR ; iPian et aD 
2006: ISoilerman et al.l20 061 with results that are consistent with 



0 10 20 30 40 50 


Rest-frame time since explosion [days] 


Fig. 7. The NIR fraction (700-2200 nm) of the quasi-bolometric flux 
(340-2200 nm) for SE SNe. (a) The values derived using optical/NIR 
data of five SE SNe are fitted separately (coloured solid lines), averaged, 
and interpolated (grey solid line and lcr contours), (b) The estimations 
from the BB fits to optical data are shown for five GRB-SNe along with 
the contours from the top panel. See main text for the references of the 
data sources. 


those for SNe 2008hw, 2009nz, and 2010ma. Moreover, the lcr 
contours (grey-shaded region) are compatible with the all BB 
estimates derived using optical photometry solely. 


The final NIR correction applied to the optical data was the 
average value between the estimates from the available NIR data 
for GRB-SNe and the estimates from the BB fits. A conserva¬ 
tive proxy of the NIR-fraction error was chosen to be the largest 
among the difference between the two estimates and their re¬ 
spective errors. Errors fluctuate between 0.07 and 0.22. We note 
that the NIR correction implies JHK magnitudes at maximum 
consistent with the upper limits presented in Figs. [Q|2] and [3] 
For instance, the brightest magnitudes derived from the NIR cor¬ 
rection are J = 22.6, H = 23.2, and K s 23.6 mag for SN 
2010ma. The corrected measurements of the quasi-bolometric 
flux are presented in Fig. [8] for the GRB-SNe 2008hw, 2009nz, 
and 2010ma. For comparison, the quasi-bolometric light curves 
(340-2200 nm) for other SE SNe are also computed and plotted 
in Fig. [8] We note that all three events lie at luminosity compa¬ 
rable to that of GRB-SNe 1998bw and 2006aj and are brighter 
than “normal” type-Ib/c SNe. Similar to the results we obtain 
for individual optical filters, SN 2010ma turns out to be brighter 
than SN 1998bw. The quasi-bolometric fluxes of SNe 2008hw, 
2009nz, and 2010ma at maximum (Table [6j are comparable to 
(1.07 + 0.07) x 10 43 erg s" 1 for SN 1998bw. 
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Fig. 8. The quasi-bolometric (340-2200 nm) light curves of SNe 
2008hw (red six-pointed stars), 2009nz (blue five-pointed stars), and 
SN 2010ma (green filled diamonds). The analytic model is shown for 
each SN with the respective colour (details in Sect. 13.31 . A handful of 
other type-Ib/c SNe are shown for comparison (see main text for the 
references). 


3.3.2. Physical parameters of the explosion 

The nickel mass M n;, the ejected mass M e j, and the ki¬ 


netic energy E k of the explosion were extracted from the lu- 



ployed the following expression to model the bolometric lumi¬ 
nosity: 


L(t) - M N j e - * 2 


Om - eco) A(c) dz + e Co 
Jo 



B(z) dz 


(1) 


where A(z) = 2ze~ 2zy+z \ B(z) = 2z e~ 2z,J+2zs+z \ x = t / Tm , y = 
T m /(2r Ni ), and .? = r m (r C o - r Ni )/(2T C oTNi)- The decay energy 
for 56 Ni and 56 Co are e N ; = 3.90 x 10 10 erg s_ 1 and ec 0 = 
6.78x10 9 erg s~| g -1 , respectively dSutherland & Wheeled 19841 
ICappellaro et aljfl997l) . The decay times are r® = 8.77 d and 
tco = 111 d. The time scale of the light curve is expressed as 


/ ^opt 

f 

(10M 3 ) 

ej 

\Pc> 


l 3F k J 


( 2 ) 


where // =* 13.8 is an integration constant ( Arne 119821) . c is the 
speed of light, and k opt is the optical opacity, which stays con¬ 
stant in time for this modelling scheme. In reality, the opacity 
depends on the composition and temperature of the ejecta, there¬ 
fore, it changes as the SN expands. Assuming a variable opacity, 
the models by iChugail (l2000i) for the bolometric light curve of 
SN 1998bw deliver an average value of 0.07 cm 2 s_j for the 
first_20 days after the explosion. The models by iMazzali et al.l 
d2000i) can reproduce the light curve of the type-Ic SN 1997ef 
at early times using a constant opacity of 0.08 cm 2 s 1 . With a 


constant opacity of 0.06 cm 2 s _1 , the synthetic light curves by 
iMaeda et al.l d2003l) manage to reproduce the data of hypernovae 
at early phases. Thus, we chose a value of k opt = 0.07 + 0.01 cm 2 
s _1 , which includes within lcr the opacity values that have been 
employed in the literature. 

Equations Q] and [2] are valid only for the photospheric phase 
(t - to £ 40 d). Given the lack of detections beyond day 40, 
no nebular c ompo nent has been considered (see appendix in 
lYalenti et al.1120081 for the complete model). The modelling 
procedure employed consists of a weighted x 2 minimisation, 
where AW, and AW / E k are free. The latter will be dubbed 
the “timescale parameter” hereafter, because it approximates the 
light-curve shape (see Sect. [4] for details). To compute M e j and 
Ek from the timescale parameter, we used the expression for the 
photospheric expansion velocity at maximum luminosity from 
[Ametl(fl982l) : 


"ph 


2Ek 

M~' 


(3) 


This quantity is critical to obtain reliable physical parameters 
of the explosion (Mazz aliet al.ll2013h . A minimum expansion 
velocity of ~ 14,000 km s -1 (SN 1998bw; IPian et alJ 20061) 
and a maximum of ~ 28,000 km s 1 (SN 2010bh; Bufano et al. 
12012 ) have been measured for GRB-SNe. Thus, we em¬ 
ployed 22,000 ± 4,000 km s 1 if estimates of the photospheric 
velocity are not available. This conservative proxy encom¬ 
passes with a 2cr confidence the photospheric velocity of most 
spectroscopically-confirmed GRB-SNe at maximum luminosity 
(see, e.g. JBufano et al.ll2012l) . 

To calculate uncertainties, we performed thousand Monte- 
Carlo simulations for each event. Assuming Gaussian errors, 
each simulation consisted of a x 2 minimisation between the 
model with a randomised opacity and the randomised quasi- 
bolometric data points. From the resulting distributions for AW, 
and M 3 /£V, we obtain the median and the standard deviation 
(lcr). Then, eq.[3]is employed to compute M c] and £ k propagat¬ 
ing the errors accordingly. When using the wide range of expan¬ 
sion velocities for SNe 2008hw and 2010ma, we computed the 
weighted average of M cj and £ k between the parameters obtain 
using the minimum and maximum photospheric velocities as de¬ 
fined above, and listed the cor respo nding lcr ranges in Table [6] 
For SN 2009nz, iBerger et al.l ( 2011 ) measure an expansion ve¬ 
locity of 17,000 km s _1 from Si n T6355, which has been identi- 
fied a s a reliable tracer of the photospheric velocity (ISauer et al.l 
120061: IValenti et al.ll2008h . Although the date of the spectrum 
(16.3 rest-frame days after the GRB) coincides quite well with 
the maximum luminosity, the spectral coverage barely extends to 
6250 A and the spectrum has low S/N. Therefore, we assigned to 
this velocity a conservative uncertainty of 1,500 km s -1 , which 
corresponds to about 30 A. The physical parameters and best-fit 
models are listed and plotted in Tableland Fig. [8] respectively. 

Figure[8]shows that the light curves are reasonably well mod¬ 
elled within the errors. In the case of SN 2009nz, however, 
IBerger et~aD (1201 lh obtain a lowe r M N i (0.35 M R of 5 6 Ni). They 
scale the /-band photometry from Cobb et all (1201 Pel) to obtain a 
V-band absolute magnitude at maximum and then they compute 
M Ni using a simplification of the formalism by Arnetti d!982l) . 
which should deliver results similar t o ours. Given that a higher 
56 Nj^mass implies higher luminosity (IColgate et al.l 19801:1Arnetti 
119821) and our data includes more flux (three GROND bands plus 
the NIR correction), our value is likely more reliable. The M e j 
quantity is consistent with that presented bv lBerger et al.l d201 ll) . 
Regarding SNe 2008hw and 2009nz, no detailed photometric 
studies have been published for these yet. 
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Table 6. Physical parameters from quasi-bolometric light curves. 


SN 

Mn i 

Mej 

E k 


^ph 

log Zy max 

^max 

Reference 


[Mo] 

[M 0 ] 

[10 51 erg] 

[10 51 M 0 erg -1 ] 

[10 3 km s~‘] 

[erg s- 1 ] 

[days] 


2008hw 

n ^ 0 + 0.08 

-0.04 

9 o +1.0 
^ -0.7 

19 + 15 

0 7+ 0 . 7 
- 0.2 

22+4 

43.1+0.3 

12 + 3 

1 

2009nz 

0.50 ± 0.04 

9 4 + 0.6 
^-0.3 

11 ±4 

i 9+0.6 

1 ‘ z '-0.2 

17 + 1.5 

43.0 + 0.2 

18 ±4 

1 

2010ma 

0 43+0-03 
-0.02 

i o +0.4 

1 -0.3 

10 + 6 

0 20 + 0 09 
-0.04 

22+4 

43.1+0.2 

10 + 2 

1 

1998bw 

0.38-0.48 
0.45 ± 0.01 

li 

3.4 ± 0.2 

50 

16 + 3 

27 

2.4 ± 0.2 

17 + 1 

43.03 + 0.01 

17.8 

2,3 

1 

2003dh 

0.25-0.45 

8 

40 

13 

~ 20 

~ 43.0 

- 18 

4 

20031w 

0.45-0.65 

13 

60 

37 

~ 18 

~ 43.2 

- 18 

5 

2006aj 

0.21 

0.26 ± 0.01 

2 

0.6 ± 0.04 

2 

3.5 ±0.6 

4.0 

0.057 ± 0.004 

19 + 1 

42.99 ± 0.02 

8.8 

6,7 

1 

2010bh 

0.21+0.03 

2.6 ± 0.2 

24 ±7 

0.7 ± 0.3 

28 

42.63 

8.0 

8,9 

2002ap 

0.10 

2.5 

4 

3.9 

14 

~ 42.4 

- 11 

10 

2003jd 

0.36 

3 

7 

3.9 

13.5 

-42.9 

- 17 

11 

2007uy 

0.30 + 0.01 

4.4 ± 0.3 

15 ± 1 

5.6 ± 1.0 

15.2 

42.83 

17.9 

12 

2008D 

0.07-0.09 

5.3 ± 1.0 

6 + 3 

25 ± 17 

~ 10 

-42.4 

- 19 

13,14 

2009bb 

0.22 + 0.06 

4.1 ± 1.9 

18 ±7 

3.8 ±5.5 

~ 20 

-42.8 

- 18 

15 


Notes. Different parts of the table correspond to GRB-SNe analysed here (top), further GRB-SNe (middle), and other SE SNe (bottom). Uncer¬ 
tainties are given at the lcr level. The u p h values correspond to measurements as defined by eq. [3] The L max and t max values are estimated at the 
maximum bolometric luminosity, where f max is defined with respect to the explosion time. 

References. (1) This pap er, (2)IIwamoto et all 1 9981 . (3 ilMazzali et al 1200 1 1 . (4)lMazzali_et_alJ|2002 , (Si lMazzali et al.l| 20()6bl (6)Kanet_alj2006 , 
(7) lMazzali etaTkOOfii (Sl lOlivares E. et all2012l ( 9)lBufano et all2012L (10 i lMazzali et alfcooi (11) 1 Valenti et all2008l (12) lRov et ali2~ , 
(13) lMazzali et alJl2008l (14) lTanaka etail2 Q09. (15) Pign ata et al.ll20Ul 


4. Discussion 

Regarding the NIR correction utilised in Sect. 3, we have to ad¬ 
dress that the extrapolation to SNe with different properties is 
the major source of uncertainty for this correction, although the 
five SE SNe selected for the analysis already cover a wide range 
of properties. A clear case of deviation from our NIR correc¬ 
tion is the contribution shown by SN 2002ap. Even though SN 
2002ap was not preceded by a GRB, it was a type-Ic event that 
showed a maximum NIR fraction of about 0.6 of the total quasi- 
bolometric flux. Moreover, the colours of SN 2010ma turned 
out to be significantly bluer before 1 < 20 d than those of SN 
1998bw. This could hint at a higher temperature of the SN en¬ 
velope and, therefore, lower NIR fluxes. Therefore, we caution 
that there might be GRB-SNe that will not fit into our estima¬ 
tion of the NIR correction. The NIR fraction could have varia¬ 
tions as large as +0.15 if we compare SN 2002ap to SN 2007uy. 
This would translate into a maximum variation of about +30% 
in the quasi-bolometric flux (equivalent to a lcr error of ~ 10%) 
and therefore in the determinations of Mm. This issue could be 
solved in the future by using a larger sample, i.e., by including 
observations of new GRB-SNe in the NIR bands. 

With the purpose of comparing the physical parameters com¬ 
puted by others for a set of different SNe, we gathered results 
from the literature in Table[6] although uncertainties were unfor¬ 
tunately not available for all. To compare the analytic method 
against the hydrodynamic al simulations, we additionally com¬ 
puted the physical parameters of the explosion for SNe 1998bw 
and 2006aj using u p h = 17,000 and 19,000 km s~' rfPian et al.l 
l2006h . respectively. In Figure [9] we plotted the kinetic energy 
per unit mass E k /M e j against the synthesised nickel mass Mm, a 
diagram that have been presented bv lBufano et al. 1 ( f201 2i). Even 
though some values have large uncertainties, we recognised a 
trend where the more energetic the SNe is, the more 56 Ni it syn¬ 
thesises (IMazzali et al.l2007l) . We note also that hydrodynamical 


(green) and analytic (blue) measurements are consistent for SN 
1998bw, despite showing significant differences for M cj and E k 
individually. This is because the ratio E k /M e j is proportional 
to u p 2 h (eq. |3}. which is a common measurement for both ap¬ 
proaches. The discrepancies in M e j and E k individually are prob¬ 
ably attributed to the different values used for the optical opacity 
and of course to the different assumptions and models employed 
(hydrodynamical or analytic). This would explain the large dis¬ 
crepancies shown for SN 2006aj as well, where the difference in 
M n ; could be explained by our inclusion of the NIR data. We 
caution that the physical parameters of the explosion might be 
highly model-dependant, especially the values obtained for M e j 
and E\-. 


5. Summary and conclusions 

Here we studied the GRB-SN connection by means of three 
individual events followed up in depth by XRT, UVOT, and 
GROND. The X-ray, UV, optical, and NIR data covered approx¬ 
imately six orders of magnitude in the radiative energy domain. 
Excluding y-ray data, this represents a very comprehensive data 
set presented for the associations GRB 081007/SN 2008hw, 
GRB 091127/SN 2009nz, and GRB 101219B/SN 2010ma. 

In Sect. |3] the light curves of the three events are thoroughly 
analysed. The host-galaxy extinction along the line of sight of 
each event is computed from the broad-band SED. The light 
curves of individual filter bands are modelled with SN 1998bw 
templates (Sect. 13. Ik The AG component is subtracted to iso¬ 
late the SN counterpart. The NIR flux was estimated from the 
data of five SE SNe and using BB fits of the optical data. This 
correction has been applied to the integrated optical flux of our 
rebrightenings to obtain quasi-bolometric light curves from 340 
to 2200 nm. We note that the NIR contribution of SN 2002ap is 
10-15% larger than that of the GRB-SNe 1998bw and 2006aj. 
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Fig. 9. The nickel mass against the envelope energy per unit mass. 
While green squares depict values obtained via hydrodynamical simu¬ 
lations, blue circles correspond to parameters measured using the ana¬ 
lytic approach explained in Sect. l3.3l or similar. When ranges are given 
in Table[6] we plotted the weighted centre of the range. We caution that 
no uncertainties are available in the literature for some measurements 
(see Table i)l. 


Moreover, the colours of SN 2010ma at early times are bluer 
than those of SN 1998bw suggesting lower NIR fluxes for this 
object. Therefore, we conclude that more NIR data is needed to 
constrain better the NIR contribution in GRB-SN light curves. 
Using an analytic model for bolometric light curves, the physical 
parameters of the SN explosion were computed for each even t 
analogous to the case of SN 2010bh in lOlivares E. et aQ d2012l) . 
We derived nickel and ejected masses of about 0.4-0.5 M e and 
1-3M 0 , respectively, and kinetic energies of about 10 52 erg, 
which are higher than those of local type-Ic SNe and compara¬ 
ble to other GRB-SN events (see Tableland Fig. [91- 

In conclusion, all three cases exhibit similarities to other 
GRB-SNe in terms of luminosity and physical parameters. SN 
2008hw turned out to be somewhat fainter and slightly bluer 
than SN 1998bw (see Table [4}. Moreover, SN 2009nz showed 
the most similarities with SN 1998bw in luminosity and evolu¬ 
tion. SN 2010ma was significantly bluer and brighter than SN 
1998bw. Both the latter and SN 2010bh have among the earli¬ 
est optical peaks ever recorded (approximately 8 days after the 
GRB) and fade more rapidly than almost every other GRB-SN, 
HN, or typical type-Ic SN. 
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Appendix A: Optical/NIR photometry 

The three Tables presented as follows are corrected for Galactic foreground extinction dSchlegel et al.|[T998b . 


Table A.l. GRB 081007/SN 2008hw. 
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Notes. The GRB trigger time is to = 54746.225 MJD. All data corrected for Ay,host = 0.68 ± 0.08 mag (SMC). Image subtraction of the host was 
performed for g'r'i'z'. (al The duration of the observation. 


Table A.2. GRB 091127/SN 2009nz. 
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1.243 

1.75 

19.23(06) 

18.99(07) 

2.079 

1.70 

20.81(04) 

20.58(04) 

20.47(01) 

20.38(02) 

2.080 

1.75 

20.09(09) 

19.86(10) 

2.198 

1.71 

20.88(03) 

20.65(04) 

20.54(03) 

20.41(03) 

2.198 

1.77 

20.09(07) 

19.90(12) 

3.207 

1.71 

21.64(07) 

21.32(04) 

21.26(05) 

21.13(05) 

3.207 

1.76 

20.54(32) 

20.52(22) 

4.212 

1.70 

22.10(09) 

21.87(06) 

21.85(10) 

21.60(07) 

4.212 

1.75 

20.88(18) 

20.99(28) 

6.175 

1.71 

22.68(13) 

22.40(10) 

22.51(14) 

22 .22(10) 

6.175 

1.75 

21.53(31) 

> 21.26 

11.10 

1.71 

23.76(12) 

22.76(06) 

22.50(10) 

22.86(13) 

11.10 

1.75 

> 21.63 

> 21.16 

23.17 

0.34 

> 22.87 

22.85(23) 

> 22.05 

> 22.15 





46.13 

3.92 

> 24.43 

24.08(13) 

23.58(27) 

23.20(22) 

46.13 

3.97 

21.74(27) 

> 21.35 

49.12 

1.70 

> 24.69 

24.49(32) 

> 23.41 

23.28(26) 

49.12 

1.75 

>21.69 

> 21.19 

54.09 

1.90 

> 21.26 

>23.88 

> 18.86 

>22.51 

54.09 

1.75 

>21.37 

> 20.93 


Notes. The GRB trigger time is to = 55162.976 MJD. Data not corrected for negligible Ayh OS t < 0.03 mag (LMC). Image subtraction of the host 
was performed for g'r'i’z'- <a) The duration of the observation. 
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Table A.3. GRB 101219B/SN 2010ma. 


t - to 

[d] 

A f a 
[ks] 

g' 

r' 

V 

z' 

t - fo 

[d] 

At 

[ks] 

J 

H 

K s 

0.342 

0.91 

19.71(05) 

19.60(05) 

19.58(06) 

19.33(07) 

0.342 

0.92 

19.10(09) 

19.02(12) 

18.73(23) 

0.356 

1.44 

19.78(05) 

19.73(03) 

19.57(03) 

19.54(04) 

0.357 

1.49 

19.26(08) 

19.08(11) 

18.84(20) 

0.376 

1.71 

19.81(06) 

19.77(02) 

19.63(03) 

19.57(03) 

0.376 

1.76 

19.33(07) 

18.98(10) 

19.07(26) 

0.446 

1.72 

20.00(06) 

19.90(04) 

19.73(03) 

19.76(04) 

0.447 

1.77 

19.61(10) 

19.31(12) 

18.74(23) 

1.353 

2.92 

21.13(09) 

21.09(04) 

20.97(05) 

20.99(07) 

1.354 

2.98 

21.06(22) 

20.64(21) 

19.88(37) 

2.381 

1.47 

21.76(17) 

21.67(07) 

21.72(09) 

21.87(15) 

6.442 

4.20 

21.44(39) 

21.06(44) 

> 19.60 

2.508 

1.61 

> 21.86 

> 22.30 

> 22.02 

> 21.83 

2.381 

1.52 

> 20.95 

21.05(30) 

> 19.11 

3.387 

1.72 

22.36(17) 

22.05(07) 

22.05(10) 

22.29(13) 

2.507 

1.52 

> 20.43 

> 19.96 

> 19.81 

5.399 

3.53 

22.64(05) 

22.26(05) 

22.27(06) 

22.53(16) 

3.388 

1.77 

>21.10 

> 20.80 

> 20.04 

6.445 

3.53 

22.99(10) 

22.49(09) 

22.46(10) 

22.78(28) 

5.399 

3.58 

>21.38 

> 20.74 

> 19.71 

9.395 

6.69 

23.17(06) 

22.41(04) 

22.23(08) 

22.93(17) 

9.377 

3.58 

>21.10 

> 20.77 

> 19.98 

10.38 

3.52 

23.07(12) 

22.49(07) 

22.34(11) 

22.83(17) 

10.38 

3.56 

>21.37 

> 21.03 

>20.21 

12.38 

3.52 

23.39(07) 

22.55(06) 

22.35(09) 

22.89(17) 

12.38 

3.57 

>21.30 

> 20.79 

>20.17 

14.39 

3.52 

23.31(17) 

22.61(15) 

22.71(17) 

22.85(28) 

14.39 

3.57 

> 20.67 

> 20.34 

> 19.13 

16.41 

6.71 

23.67(19) 

22.83(05) 

22.46(16) 

22.85(24) 

16.40 

3.59 

> 20.95 

> 20.38 


24.37 

4.18 

> 23.72 

23.23(15) 

22.72(18) 

23.00(20) 

24.37 

4.23 

>21.28 

> 20.67 

> 19.35 

26.35 

0.78 

> 22.36 

> 22.79 

> 21.57 

> 21.59 

26.35 

0.83 

> 20.35 

> 20.11 

> 19.50 

26.43 

9.85 

> 23.95 

23.72(16) 

22.77(13) 

22.84(19) 

26.44 

4.94 

>21.71 

> 21.19 

> 19.84 

35.44 

4.88 

> 24.46 

23.75(11) 

23.28(14) 

23.60(34) 

35.44 

4.93 

>21.44 

> 21.08 

> 19.47 


Notes. The GRB trigger time is t 0 = 55549.686 MJD. All data corrected for Av- ihost = 0.12 ± 0.01 mag (SMC). (a> The duration of the observation. 


Appendix B: Sequences of standard stars 

The sequence of reference stars in the held of GRB 091127/SN 2009nz are taken from lFilgas et alj d20 1 ll) . Stars from the 2MASS 
catalogue JSkrutskie et al.H2006t) are used for the JHK S bands. 


Table B.l. Reference stars in the field of GRB 081007/SN 2008hw. 


RA 

[°] 

Dec. 

[°] 

g' 

r' 

V 

z' 

339.91132 

-40.14995 

19.334 + 0.013 

18.940 + 0.011 

18.795 + 0.016 

18.658 + 0.017 

339.91385 

-40.15378 

19.050 + 0.011 

18.495 + 0.007 

18.274 + 0.012 

18.075 + 0.012 

339.92912 

-40.17095 

14.897 + 0.001 

14.399 ± 0.001 

14.220 + 0.001 

14.091 +0.001 

339.94805 

-40.15443 


19.887 ± 0.020 

18.788 + 0.016 

18.285 + 0.014 

339.95393 

-40.11101 

19.971 +0.020 

19.147 + 0.011 

18.807 + 0.016 

18.641+0.018 

339.97304 

-40.11895 

19.890 + 0.019 

19.201+0.012 



339.97329 

-40.17987 

19.631+0.016 

18.129 + 0.006 

16.949 ± 0.005 

16.399 ± 0.004 

339.98310 

-40.12094 

19.542 + 0.015 

18.780 + 0.009 

18.488 + 0.013 

18.309 + 0.014 

339.99130 

-40.17932 

17.947 + 0.005 

17.691 +0.004 

17.588 + 0.007 

17.517 + 0.007 

Notes. All are observed magnitudes 

in the AB system. 





Table B.2. Reference stars in the field of GRB 101219B/SN 2010ma. 

RA 

[°] 

Dec. 

[°] 

g' 

r' 

V 

z' 

12.22124 

-34.52946 

22.85 + 0.10 

19.87 ± 0.05 

18.365 + 0.019 

17.720 + 0.018 

12.25499 

-34.54326 

19.021 +0.031 

18.006 + 0.011 

17.568 + 0.012 

17.368 + 0.015 

12.26985 

-34.56368 

21.02 + 0.07 

19.519 + 0.034 

18.475 + 0.022 

18.088 + 0.025 

12.22871 

-34.56753 

15.350 ± 0.003 

14.928 ± 0.002 

14.772 ± 0.002 

14.674 ± 0.003 

12.24608 

-34.57123 

17.328 + 0.010 

17.099 ± 0.006 

17.017 + 0.009 

16.970 + 0.011 

12.21783 

-34.57419 

20.41 ± 0.06 

19.226 ± 0.029 

18.663 + 0.026 

18.484 + 0.039 

12.26934 

-34.58304 

16.972 ± 0.007 

16.599 ± 0.004 

16.446 ± 0.006 

16.398 ± 0.007 

12.25056 

-34.58787 

16.996 ± 0.007 

15.583 ± 0.003 

14.623 ± 0.002 

14.166 + 0.002 

12.21137 

-34.59066 

18.380 + 0.017 

18.173 + 0.013 

18.097 + 0.016 

18.125 + 0.027 


Notes. All are observed magnitudes in the AB system. 
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Appendix C: Blackbody fits 

Here we present the blackbody fits for the analysed GRB-SNe, which were utilised to estimate the NIR contribution (Sect. [33711 . 






Fig. C.l. Blackbody fits to the optical photometry of SN 2008hw. Colour temperatures are about 5,000 K. Points with only a lower error bar are 
upper limits. The blue shaded region shows the area between the lcr contours, where the central line is the best fit. 






Fig. C.2. Blackbody fits to the optical photometry of SN 2009nz. Colour temperatures evolve from ~ 7, 000 to ~ 4, 000 K approximately. Point, 
line, and region coding are the same as in Fig. 1C. II 
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Fig. C.3. Blackbody fits to the optical photometry of SN 2010ma. Colour temperatures evolve from ~ 6,000 to ~ 4,000 K. Point, line, and 
region coding are the same as in Fig. 1C. II 










































